Origin of four-fold anisotropy in square lattices of circular ferromagnetic dots 
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We discuss the four- fold anisotropy of in-plane ferromagnetic resonance (FMR) field H r , found 
in a square lattice of circular Permalloy dots when the interdot distance a gets comparable to the 
dot diameter d. The minimum H r , along the lattice (11) axes, and the maximum, along the (10) 
axes, differ by ~ 50 Oe at a/d = 1.1. This anisotropy, not expected in uniformly magnetized 
dots, is explained by a non-uniform magnetization m(r) in a dot in response to dipolar forces in the 
patterned magnetic structure. It is well described by an iterative solution of a continuous variational 
procedure. 

PACS numbers: 75.10.Hk; 75.30.Gw; 75.70.Cn; 76.50.+g 
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Magnetic nanostructures are of increasing interest for 
technological applications, such as patterned recording 
media [lj , or magnetic random access memories 2] . One 
of the most important issues for understanding their col- 
lective behavior is the effect of long-range dipolar inter- 
actions between the dots [3J . For the single-domain mag- 
netic state of a dot, the simplest approximation is that 
dots are uniformly magnetized and interactions only de- 
fine relative orientation of their magnetic moments |j] . If 
so, the system of dipolar coupled dots in a square lattice 
should be magnetically isotropic. 

However, in all known experimental studies of closely 
packed arrays of circular dots, a four-fold anisotropy 
(FFA) was found, either by Brillouin light scattering 
[5j , ferromagnetic resonance (FMR) [fjfor magnetization 
measurements (from hysteresis loops) ]3, ■ It is impor- 
tant to note that FFA exists in both unsaturated samples 
and saturated ones (i.e. above vortex annihilation point 
on the hysteresis loop). Hence it cannot be only associ- 
ated with vortex formation suggested in Ref. |2] . It was 
instead qualitatively related to stray fields from unsatu- 
rated parts of magnetization inside the dots [f| . However 
no quantitative description of FFA in such systems was 
given up to now. So the aim of this study is to explain 
quantitatively the deviations from isotropy in terms of 
modified demagnetizing effect in a patterned planar sys- 
tem at decreasing inter-dot distance, from the limit of 
isolated dot to that of continuous film. The choice of X- 
band FMR techniques for this study has an advantage in 
eliminating possible interference from domain (vortex) 
structure |9j. The variational theoretical analysis is fol- 
lowed by micromagnetic simulations. 

Permalloy (Py) dots were fabricated with electron 
beam lithography and lift-off techniques, as explained 
elsewhere 101. The dots of thickness t = 50 nm and 



diameter d — 1 /mi were arranged into square arrays 
with the lattice parameter a (center to center distance) 
varying from 1.1 /im to 2.5 /im. The dimensions were 
confirmed by atomic force microscopy and scanning elec- 
tron microscopy. Room temperature FMR studies were 
performed at 9.8 GHz using a standard X-band spec- 
trometer. The dependence of the FMR field H r on the 
azimuthal angle <fn of applied field H with respect to the 
lattice [10] axis for almost uncoupled dots (a = 2.5 fan) 
is shown in Fig. Only a weak uniaxial anisotropy 

of H r ((pii) is present here, which can be fitted by the 
simple formula H r (ifH) = H r ^ av + H2Cos2fH- For the 
a = 2.5 /im sample, we found the average peak posi- 
tion H r av ks 1.13 kOe and the uniaxial anisotropy field 
H.2 ~ 5 Oe. The latter value remains the same for the 
rest of our samples, so this uniaxial anisotropy is most 
probably caused by some technological factors. 

With decreasing distance a between dots, two changes 
are observed in the -ff r (</?#) dependence. First, H r ^ av 
decreases to « 1.09 kOe at a = 1.1 /im (Fig. Wp)- Second, 
a four-fold anisotropy (FFA) is detected in the samples 
with a <1.5 /im by pronounced minima of H r (ifH) at fn 
close to the lattice (11) axes. This behavior is fitted by 
H r (a,(p H ) — Hr iav (a) + H 4 (a) cos4ip H + H 2 cos2(^ H , as 
shown in Fig. The interdot distance dependence of 
H r ^ av and FFA field H4 is shown in Fig. Also such 
anisotropy is detected in the FMR linewidth, smaller for 
(11) than for (10) case (reaching ~30% at a/d = 1.1). 

The FFA effect, which could not arise in uniformly 
in-plane magnetized cylindrical dots, is evidently re- 
lated to a non-uniform distribution of the magnetization 
m(r, <f, z) (in cylindric coordinates < r < R = d/2, < 
if < 2-7T, < z < t). A similar effect was discovered using 
Brillouin light scattering and magnetization reversal 
in such systems under weak enough external fields, 
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FIG. 1: In-plane FMR field in square lattices of 1 /im circular 
Py dots as a function of field angle tpn- a) The data for lattice 
parameter a — 2.5/im are well fitted by uniaxial anisotropy 
(solid line), b) At a = l.ljum, the best fit (solid line) is 
a superposition of FFA and uniaxial anisotropy (separately 
shown by dashed line). 
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FIG. 2: a) Average FMR field H r ,av and FFA field H 4 as func- 
tions of interdot spacing a. The points are the experimental 
data and the solid lines present the 1st iteration theory (the 
limits mark H r of isolated dot and continuous film) . b) The 
same data plotted against (d/a) 2 for H r ^ a v and (d/a) 4 for H4 
give excellent linear fits (dashed lines). 



(see Appendix) 

£ = -H.m + ^£^M|m g .g| 2 , (1) 

a 9 

where f(u) = 1 — (1 — e~") /u Q and the vectors 
of the 2D reciprocal lattice are (^-dependent: g = 
(27r/a)(ni cos ifiH — n-i sunpn, n i shi(p# + ri2 cos ipn) (for 
H || x and integer ^1,2)- The variation of exchange 
energy at deformations on the scale of whole sample 
is of the order of stiffness constant (~ 10~ 6 erg/cm 
for Py) and it can be neglected beside the variation 
~ HM s d 2 ~ 10~ 2 erg/cm of terms included in Eq. ^ 
If the dot magnetization has constant absolute value: 
m(r) = M s (cos <p(r), sin <p(r)), its variation: Sm(r) = 
z x m(r)6ip(r) (where z is unit vector normal to plane), is 
only due to the angle variation 6<p(r). Using the Fourier 
transform <5m g = z x ^ , m s _ s >5ip s > in the condition 
6E = leads to the equilibrium equation for the Fourier 
amplitudes: 

(m g _ g , x g') • z. (2) 
It can be suitably solved by iterations: 

(„) _ **_ V- f(9't) ( (n-l) t) 

(m^xgO-z, (3) 

starting from uniformly magnetized dots as zeroth it- 
eration: mg°j = 0, TOg°z = 2ttRM s Ji(gR)/g (with the 
Bessel function Ji). Already the 1st iteration (including 
the inverse Fourier transform): 



which displace vortices in each dot. This can be modeled 
by displacements of two oppositely in-plane magnetized 
uniform domains [ll) . But in the presence of external 
fields strong enough to observe FMR, one has to assume 
a continuous (and mostly slight) deformation of m(r, tp). 
The simplest model for such deformation uses a varia- 
tional procedure with respect to a single parameter [l2j| . 
However, as will be shown below, the non-uniform mag- 
netic ground state of this coupled periodic system results 
from a rather complicate interplay between intra-dot and 
inter-dot dipolar forces, which requires a more general 
variational procedure. 

Assuming fully planar and z-independent dot mag- 
netization with the 2D Fourier amplitudes m g = 
J e lg r m(r)(ir, the total (Zeeman plus dipolar) magnetic 
energy (per unit thickness of a dot) can be written as 



Ji{gR) cos(g-r), (4) 

(with the Heavyside 9 function) reveals the FFA be- 
havior, due to the rotationally non- invariant product 
9x9y The calculated maximum variation of ip(r) — 
arcsin[m y (r)/M s ] in (10) field geometry is ~ 60% big- 
ger than in the (11) geometry (Fig. |3 upper row). This 
expected behavior persists upon further iterations. Our 
analytic approach was checked, using the micromagnetic 
OOMMF code 0| ° n a 9x9 array of considered disks 
(Fig. lower row) at standard values of M s = 0.83 kOe 
and exchange stiffness 1.3 • 10 -6 erg/cm for Py. The 
distributions obtained in this way for the central disk in 
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FIG. 3: Density plots of the equilibrium magnetization 
angle <p(r) for two field geometries, ipH = ((10)) and 
tpH = tt/4 ((11)). Upper row: calculated from the sum, Eq. 
0] over 100x100 sites of reciprocal lattice at parameter values 
a = 1.1 /im, H — 1.1 kOe, M s — 0.83 kOe. Lower row: mi- 
cromagnetic calculation by OOMMF code for the central disc 
in the 9x9 array. 



the array are within ~ 10% to the analytic results of the 
1st iteration. 

The FMR precession of m(r) is defined by the internal 
field Hj(r) = H + h(r) through the local dipolar field 



h a {r) 



4-7T 



— 2^ 2 m g ,/3 cos (g • r) , (5) 



N XlV (r) dr 



r<R 



= {^/a 2 )Y,fWl{gR)/g\ (7) 

and N z = 4tt — 2N X . At a — > oo, they tend to the single 
dot values [HI which are for t/R = 0.1: « 0.776 

and N ( z d) w 11.01. Using N^ d) instead of A^(r) in 
Eq. El accurately reproduces the single dot FMR limit 
Hr ~ 1.15 kOe (estimated from Fig. [2h>) - Otherwise, 
for decreasing interdot distance, a — > d, the 1st itera- 
tion values, Eq. used in Eq. H3 well describe the 
tendency of H r _ av (a) towards the continuous film limit 
HP = y/H$ + 4tt 2 M s 2 - 2ttM s w 0.96 kOe (Fig. Hi). 

Finally, by calculating the true local FMR fields H r (r) 
from Eqs. El and El the field dependent absorption is 
obtained as 1(H) oc f r<R 5(H — H r (r))dr. Then the 
FMR fields, H r defined from maximum of 1(H) in two 
geometries, display FFA in a good agreement with the 
experimental data (Fig. 0). This effect is due to the 
fact that stronger deformation of magnetization stronger 
suppresses the demagnetizing effect (the differences N z — 
N x ^ y ) and thus enhances H r . Also it produces a bigger 
spread of local resonance fields H r (r) and thus broadens 
the FMR line, again in agreement with our observations. 

In conclusion, it is shown that under in-plane magnetic 
fields, H, even strong enough for FMR, the dipolar cou- 
pling in a dense lattice of circular magnetic dots is able 
to produce a continuous deformation of the dot magne- 
tization, strongest for the field orientation along lattice 
axes. 
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(a, j3 = x,y). The 1st iteration for h(r) corre- 
sponds to the zeroth iteration for m(r) = (M s , fi y , fx z ), 
which now includes the uniform FMR amplitudes 
fiy, [i z . Then the local demagnetizing factors N x (r) — 
-h x (r)/M s , N VtZ (r) = -h VtZ (r)/ n VtZ define the local 
FMR field H r (r): 

H r (r) = \jHl + M*[N z (r)-N y (r)Y/A 

- M s [N z (v)+N y (v)-2N x (v)]/2 (6) 

(here H rj 3.3 kOe). The average FMR field is defined 
by the isotropic averaged demagnetizing factors 
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APPENDIX 



where the Fourier amplitudes of the dipolar field are: 
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dr / dz'e i(s r+<?z,) / dr' / dz" 

J-oo J J-t/2 

d 2 m p{ v ') 



t/2 



E 



dr a drp 



13 — y / |r-r'| 2 + (z'-z") 2 
To calculate them, we express the lattice magnetization 
m^(r') through its Fourier amplitudes: 



mp(r') 



i 2 ^ 
s' 



and then introduce the factor e^ s ' r ~i z ) into the integral 
in dr' dz' , and the compensating factor e~ l ( g r -<? z ) into 
the integral in drdz" . Then the spatial integrations in 
Ed are done accordingly to the formulas: 



For fully planar and z-independent dot magnetization, 
the dipolar energy per unit thickness of a dot in the lat- 
tice is: 



E d = — dz dz' dr dr'y^m a (r) 

Zt J -t/2 J -t/2 Jc J 

mp(r') 



a.,/3 



d 2 

dr a dr/3 



^\r-r'\ 2 + (z-z'f 



where the 2D integrations J dr and J dr are respectively 
over the unit cell and over the entire plane. It can be 
also presented as 



d re i(g-g')r = a 2 S g ^, 



*/2 ft/2 , . 

dz dzV^"- 2 ) = 4 sin -^ 

-t/2 J-t/2 Q 2 



dz' / dr'e iB - ( - r - r " )+lq{z '- z " ) 
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'dr*dr f3 J lr _ r , l 2 + {z ,_ z , f) 2 



g 2 + q 2 



Finally, the momentum integration 
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sin 2 (gi/2) 
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leads to the result included in Eq. 



